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ABSTRACT 

Recent High Energy Stereoscopic System (HESS) observations show that microquasars in high- 
mass systems are sources of very high energy 7-rays. A leptonic jet model for microquasar 7-ray 
emission is developed. Using the head-on approximation for the Compton cross section and taking 
into account angular effects from the star’s orbital motion, we derive expressions to calculate the 
spectrum of 7 rays when nonthermal jet electrons Compton-scatter photons of the stellar radiation 
field. The spectrum of Compton-scattered accretion-disk radiation is also derived by approximating 
the accretion disk as a point source of radiation located behind the jet. Numerical results are compared 
with simpler expressions obtained using d-function approximations for the cross sections, from which 
beaming factors are derived. Calculations are presented for power-law distributions of nonthermal 
electrons that are assumed to be isotropically distributed in the comoving jet frame, and applied to 
7-ray observations of LS 5039. We conclude that (1) the TeV emission measured with HESS cannot 
result only from Compton-scattered stellar radiation (CSSR), but could be synchrotron self-Compton 
(SSC) emission or a combination of CSSR and SSC; (2) fitting both the HESS data and the EGRET 
data claimed to be associated with LS 5039 requires a very improbable leptonic model with a very 
hard electron energy distribution. Because the 7 rays would be variable in a leptonic jet model, the 
data sets are unlikely to be representative of a simultaneously measured 7-ray spectrum. We therefore 
attribute EGRET 7 rays primarily to CSSR emission, and HESS 7 rays to SSC emission. Detection 
of periodic modulation of the TeV emission from LS 5039 would favor a leptonic SSC or cascade 
hadron origin of the emission in the inner jet, whereas stochastic variability alone would support a 
more extended leptonic model. The puzzle of the EGRET 7 rays from LS 5039 will be quickly solved 
with GLAST. 

Subject headings: gamma rays: microquasars—radiation processes: nonthermal 


1. INTRODUCTION 

X-ray binaries with jets, or microquasars, are com¬ 
mon in our Galaxy, wi t h sa 16 now known (for a re¬ 
cent review see, e.g.. lParedegll2005il . About one-third are 
high-mass X-ray binaries (HMXBs), including Cygnus 
X-l, Cygnus X-3, LS 5039, and LSI+61°303, and the re¬ 
mainder are low-mass X-ray binaries (LMXBs), including 
GRS 1915+105, GRO J1655-40, Sco X-l, and IE 1740.7- 
2942. The compact companions are a mixture of black 
holes and neutron stars, and the radio activity of the mi¬ 
croquasars is about equally divided into persistent and 
transient behaviors. 

Recent observations llAharonian et al ] [20f)l made 
with the High Energy Stereoscopic System (HESS) show 
that the high-mass microquasar LS 5039 is a source 
of very high energy (VHE) 7 rays in the ~ 200 GeV 
- 10 TeV range, confirming its earlier tentative iden¬ 
tification with the EGRET source 3EG J1824-1514 
( Paredes et al ) '200 0:. A second high-mass microquasar 
system, LSI+61°303 (V615 Cas) , is associated with 
the COS-B source 2CG 135+01 l|Hermsen et all Il977t 
iGregorv fc Tavloi l I1978D and the EGRET source 3EG 
J0241+6103 (Kni ffen et alJ 1199^1 . but is too far north 
for observations with HESS. The EGRET source 3EG 
J1824-1514 associated with LS 5039 shows marginal evi¬ 
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dence for variability (iTorres et al.l2003il . In contrast, the 
EGRET light curve of 3EG J0241+6 103, the counterpar t 
to LSI+61°303, is stron gly variable llTavani et al .11 19981) . 
Moreover, iMassil 12004H performed a timing analysis of 
the 3EG J0241+6103 data and found a most probable 
period of 27.4 ± 7.2 days, compared to its 26.5 day or¬ 
bital period. 

Evidence for stochastic and periodic variability of these 
sources at 7-ray energies would argue in favor of a lep¬ 
tonic microquasar jet model similar to blazar jet models 
(for r ecent reviews of microquasa r models, see iRomerol 
2005; Fender fe Maccaronel li2004L especially if the X- 
ray an d 7-ray em issions display correlated variability 
iGupta et al.ll2005l L The importance of Compton scat¬ 
tering of external photons to prod uce gamma-rays in 
blazar jets was first considered by iBegelman fc S ikora 
ifTosl and iMelia fc Konigll 1198911 and la ter, in view of 
the C omp ton Observatory dis coveries, by iDermer et all 
fl99l and lSikora et all (|1994f) . A microquasar jet model 
differs importantly from a blazar jet model through 
the addition of the stellar radiation held from the 
high-mass star and the periodic orbi tal mod ulat ion of 
the binary stellar system iGeorganopoulos et a i m , 
fc^man Bernado et al . 200211 . Although iCasares et all 
12005ft claim low significance periodic variability when 
folding the HESS data for LS 5039 with its orbital pe¬ 
riod, the EGRET data showed no comp elling evidence 
for either stochastic or periodic variability iParedes et all 
HU- The X-rays from LS 5039 are, however, moder¬ 
ately variable. RXTE observations in the 3-30 keV 
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range may show periodic variability correlated with the 
periastron passage of LS 5039, so it is unceratin whether 
the X-rays are associated with accretion disk or the jet 
dBosch-Ramon et al.feooli l. 

In addition to stellar and accretion-disk emis¬ 
sions, microquasar emission from the jet will pro¬ 
duce a variable multiwavelength continuum consist- 
ing of radio/IR and jet X-ray ( Markoff et al . 1200ll) 
synchrotron radiation. Nonthermal 7 -ray emis¬ 
sion is likely to_ori ginate from sy nchrotron self- 

Compton (SSC) (lAtovan fc Aharonianl Il999t) and ex- 
ternal Compton (EC) p roces ses dLevinson fo Blandfordl 
Il996t iCeorganopoulos et al.112002) by these same jet 
electrons. The bright high-mass star makes an im¬ 
portant contribution to the external radiation field in 
HMXBs, whereas the accretion di sk is the dominant ex - 
ternal photon source in LMXBs llCrenier et al.ll2C)05a ). 
Compton-scattering leptonic jet models of the 7 -ray 
e missio n from LS 5039 _and LSI+ 61°303 are presented 
bv iBosch-Ramon fe Paredesi (l2004albD . 

In this paper, we perform a Compton-scattering anal¬ 
ysis of the jet 7 -ray emission from HMXB microquasars 
for a leptonic jet model, focusing on Compton-scattering 
effects from the azimuthal variations of the stellar radi- 
ation field using parameters inferred from observations 
llMcSwain et alJl200H ICasares et alJ[2005 1 of LS 5039, 
which has a period of 3.90603T0.00017 days. We assume 
that the twin jets of the microquasar are oriented nor¬ 
mal to the orbital plane of the compact object and star; 
geometrical complications of precessing jets are not con¬ 
sidered here. The orbital variations of the bright O or B 
stars introduce interesting kinematic variations that ap¬ 
pear in the 7 -ray emission spectrum if the 7 rays are due 
to stellar photons that are Compton-scattered by non¬ 
thermal jet electrons, including variations of peak vF v 
photon energy and inferences of the locations of the 7 - 
ray emission site. This emissi on is also subject to th e ef¬ 
fects of 77 absorption llBottcher & Dermeril2005t IDubusl 
120051) . although this effect is not included in the calcula¬ 
tions shown here. 

Angle-dependent effects on Compton-scattered jet ra¬ 
diation are treated in Section 2. Approximations made in 
the derivation are clearly enumerated, so that they can 
be relaxed in more detailed numerical treatments. In 
particular, the analysis employs a fixed electron distri¬ 
bution. Spectral calculations are presented in Section 3 
using parameters appropriate to LS 5039. Difficulties to 
fit the combined EGRET and HESS spectra of LS 5039, 
if assumed to be simultaneously radiated, are discussed 
in Section 4. Implications for establishing the nature of 
microquasar 7 -ray emission from LS 5039 from further 
HESS and upcoming GLAST observations are also con¬ 
sidered. A summary of the results is given in Section 
5. 

2. GAMMA RAYS FROM COMPTON-SCATTERED 
STELLAR RADIATION 

The geometry of the microquasar system is shown in 
Fig. 1. The generic system considered here is a HMXB 
with parameters taken from observations of LS 5039, but 
our results are also applicable to LMXBs when the accre¬ 
tion disk is approximated by a hot spot at the base of the 
twin jets. The star and the compact object, separated by 
distance d , are assumed to follow circular orbits around 




Fig. 1.— Geometry of scattering in the microquasar system, 
(a) Stationary frame K* of the microquasar system. The stellar 
azimuthal angle </>* = 0 when the high-mass star is closest to the 
observer. The photon azimuth in the comoving frame is then </>* — 
7 r. (b) Comoving frame K' of the jet plasma. In the head-on 
approximation, the photon is scattered along the direction of the 
scattering electron. See text for definition of symbols. 


their common center of mass. Material from Roche-lobe 
overflow in low-mass systems, or stellar winds in high- 
mass systems, forms an accretion disk surrounding the 
compact object. As the matter accretes onto the neu¬ 
tron star or black hole, plasma is assumed to be ejected 
transverse to the orbital plane in the form of twin jets. 

2 . 1 . Model Assumptions 

Given that the dominant radiation mechanism for 
microquasar 7 -ray emission is undecided, for example, 
wheth er microquasar emiss i ons have a primary hadr onic 
(e.g., iR.omero et all 120031 lAharonian et all l20Q5aD or 
lepto nic (e.g.. IBosch-Ramo n et al.l l2005b: Pa redes et alJ 
120051) origin, an overly complicated leptonic jet model 
for microquasar emission seems unjustified at this time. 
A number of simplifying assumptions are therefore made 
that can be relaxed in more detailed treatments. Specif¬ 
ically, we assume that 

1 . the jets are steady and oriented in the direction 
normal to the orbital plane of the binary system, 
with the jet electrons confined to the axis defining 
the jet direction; 

2 . the orbit is circular; 

3. the star and accretion disk can be treated as point 
sources of radiation; 

4. the nonthermal electrons are isotropically dis¬ 
tributed in the comoving frame with a fixed energy 
distribution, and most of the radiation is emitted 
at a specified distance above the orbital plane; and 

5. cascade processes can be neglected. 

Assumption 1 means that we neglect additional peri¬ 
odicities of the jet ejection process possibly associated 
with periastron passage, or disk precession that could 
lead to additional modul ation of the q - ray si gnal. This 
latter periodicity (e.g., iKaufman Bernado et al.l 120021 
iTorres et alJl2005l) could be different from the orbital pe¬ 
riodicity. Moreover, the jets need not be aligned normal 
to the orbital plane. For example, the jets in V4641 Sgr 
may lie within ~ 36° of the orbital plane llButt et alJ 
12003 1. These effects can be dealt with in the present 
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model by treating a more complicated geometry, and X- 
ray and TeV observations should be cross-correlated to 
search for periods unrelated to orbital motions. Helical 
motions of jet traject ories, which may operate in black- 
hole jet systems (e.g., lRiegeill2004D . would also require a 
more detailed treatment than considered here. 

Regarding assumption 2, scattering kinematics for 
highly eccentric orbits can be developed on the basis 
of the treatment presented here. LS 5039 executes a 
moderately elliptical orbit in a system with a 22.9 Mq 
star of spectral ty pe 06.5, with eccentricity = 0.35 
(jCasares et, all 1200511 . This ef fect will be displaye d in 
both the scattering kinemat ics dParedes et, alJl2005H and 
77 absorption calculations llDubusll2005l) . but does not, 
to first order, change the conclusions we draw by study¬ 
ing anisotropic Compton scattering and 77 absorption 
for circular orbits. 

For the purposes of modeling VHE 7 -rays, analysis 
of a microquasar system must take into account the 
angle-dependence of the stellar and accretion-disk radia¬ 
tion field both in the Thomson and Klein-Nishina (KN) 
regime s of scattering. Although iGeorganopoulos et alJ 
( 2001 1 provide a useful method to calculate scattered 
radiation spectra in the KN regime, it is applicable as 
presented only to surrounding isotropic radiation fields, 
and for jets w i th fix ed iso tropic electron distr i butions . 
iParedes et~all (l2005|) and lBosch-R,amon et alJ l|2005blf 
also treat the anisotropy of the stellar radiation field in a 
simplified fashion, without fully taking into account the 
effects of the directional target photon field. 

In this work, we treat the angular dependence of 
Compton-scattered stellar and accretion-disk radiation 
fields employing an accurate approximation to the Comp¬ 
ton cross section in the head-on approximation that 
is valid throughout the Thomson and KN regimes. 
To avoid over-complicated expressions, the star and 
accretion-disk are treated as point sources of photons 
scattered by nonthermal electrons in a relativistically- 
moving jet with negligible extent (assumption 3). These 
assumptions simplify the derivation, but in no way vi¬ 
olate the essential geometry of the microquasar system, 
and are straightforward (if tedious) to relax. Neverthe¬ 
less, it is important to note that the companion star in 
a HMXB can subtend as much as ~ 10% of the full sky 
as seen from a location close to the compact object. 

We also assume a non-evolving electron distribution 
and calculate emission from a fixed jet location (assump¬ 
tion 4). As the jet plasma moves away from the micro¬ 
quasar system, the electrons are energized, for example, 
through internal collisions of ejected plasma shells or via 
shocks formed during interactions of the jetted plasma 
with the external medium. The nonthermal electrons 
subsequently lose energy through expansion and radia¬ 
tive losses. In more detailed treatments, electron-energy 
evolution through synchrotron, Compton, a nd adi abat ic 
losses are considered (e.g. , lAtovan fc Aharonianl 119991: 
iBosch-Ramon et al.l l2005at iGupta et al.l 120051). When 
stellar radiation fields are important, a consistent treat¬ 
ment of the evolution of the nonthermal electron Lorentz 
factor distribution with location requires, in addition to 
the synchrotron, SSC, and adiabatic energy-loss rates, 
electron energy-loss rates obtained from the Compton- 
scattered stellar and accretion-disk radiation fields, as 
derived here. Given that 7 -ray telescopes integrate over 


a long time to accumulate signal in comparison with 
Compton and synchrotron cooling time scales for the 
7 -ray emitting electrons, as is easily demonstrated, the 
detected emission is well-approximated as the radiation 
from a time-averaged electron distribution. Yet the 
power-law form of the electron distribution remains an 
assumption that may not agree with calculations of elec¬ 
tron evolution. 

Finally, we do not treat the electromagnetic cas¬ 
cades formed through 77 attenuation as 7 -rays propa¬ 
gate through the anisotropic radiation field of the jet 
an d sta r (see, e.g ., P rot heroe .eJLalJ llM^ tBed na re k jjllal 

lAharonian et al.ll2005at) . Cascading effects could be im¬ 

portant when the absorption depth to 77 pair produc- 
tio n attenuation is large, n amely within the inner jet 3 
llBottcher fc Dermei! l2005t iDubusI 12005 1. This effect 
could, in principle, reduce the level of modulation of 
the VHE 7 -ray signal from 77 attenuation and introduce 
new spectral components. Relaxing assumption (5) in¬ 
troduces, however, new assumptions about the strength 
and geometry of the magnetic field in the vicinity of the 
microquasar outside the jet plasma that are difficult to 
constrain. 


2 .2. Stellar Radiation Field 

We treat a system (Fig. la) where the jet outflow 
has constant bulk Lorentz factor T = 1/ ^/l — (3 2 . In 
the stationary frame of the microquasar system (where 
the center-of-mass of the compact object and high-mass 
star is at rest), the plasma will have traveled a distance 
x = pet * from the compact object and reached a dis¬ 
tance r = \Jx 2 + d? from the companion star after time 
t*, measured from the moment of ejection of the jetted 
plasma at x = 0 . 

Let (t>*(t) = w»f» represent the star’s orbital phase 
measured from the angle when the high-mass star is 
closest to the observer, where w* = 2tt/P is the star’s 
angular frequency and P is the period (from hours to 
days for LMXBs, and days to weeks for HMXBs). Pho¬ 
tons from the companion star impinge on the outflowing 
jet plasma at an angle d* = arccos/I* to the jet axis 
as measured in the stationary frame of the star, where 
p* = x/r = 1/ yj\ + d 2 /x 2 . The angle 9 = arccos \i is the 
inclination of the observer with respect to the jet axis. 

The energy flux at distance r from a uniform brightness 
sphere (i.e., the star) of radius R* is 

d£* _ 7 tRI 

dAdtde * r 2 e * ’ 


where £* is the radiated photon energy and e* = 
hv*/m e c 2 is the dimensionless photon energy. The in¬ 
tensity 


B 


/“(©) 


2 ro e c 3 e 3 

Ac[exp(e*/0) - 1] 


for a blackbody, where Ac = h/m e c = 2.42 x 10 ~ 10 
cm is the Compton wavelength of the electron and 0 = 
fcs?*/TO e c 2 is the dimensionless temperature of the star. 

Consider a star with luminosity L* and temperature 
T*, so that the stellar radius R* = ^Lx/AttctsbT? , and 


3 Inner jet refers to jet locations x < d. 
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cxgB is the Stefan-Boltzmann constant. The differential 
energy density 


<b(e *;0 


1 d£* 

c dAdtde * 


Hence 


l bb 


(e*,H*;r) = ■ 


- y *)<$(</>* - </>♦) 

exp(e*/ 0 ) — 1 


(1) 


where 

o 15L* 

7 / = - 

* 47 r 5 c 0 4 r 2 

The stationary-frame photon density n* h (e*,n*) = 

u*(e*, fi*)/(m e c 2 e*). 

The f* = i/E* spectrum of the star measured by an 
observer located a (luminosity) distance away from 
the star is 


= 15 U e 4 

/e 47 r 5 0 4 d 2 exp(e/0) - 1 ’ 


( 2 ) 


where e = e*/(l + z) = e* for the low redshift (z <C 1 ) 
sources considered here. 


2.3. Stellar Radiation Photons Compton-Scattered by 
Jet Electrons 


The detailed derivation of the Compton-scattered stel¬ 
lar radiation (CSSR) spectrum is given in Appendix 
A. The vF v spectrum resulting from Compton-scattered 
stellar radiation (CSSR) for a uniform emitting region 
filled with an isotropic comoving distribution of electrons 
is 


e C* 


3co- T M^D e2 

32-7T d 2 



K( 7) 


x 


« 2 [9 2 (y + 7)^1 - 


1. 

y' 


2e 


S D ^ybg 


h+ &C (3) 


where y = 1 - (e/< 5 D 7 ), 


d D = [T(l-/3 M )]- 1 , (4) 

is the Doppler factor, e = ho/m e c 2 is the dimensionless 
observed photon energy, TV'( 7 ) is differential distribution 
of electrons with comoving Lorentz factors 7 , and d^ is 
the luminosity distance to the source. The functions a, 
b, 1 1 , I 2 and I 3 are defined in Appendix A. 


2.4. Compton Spectrum from Point Source of Radiation 
Field Behind Jet 

The accretion disk provides a source of external ra¬ 
diation that enters the jet from behind. For a Shakura- 
Sunyaev accretion disk, the distance where the transverse 
extent of the accretion disk can be neglected and the ac¬ 
cretion disk can be approximated by a point source is 
given by 

x » T 4 r ff (5) 

llDermer fo SchlickeiseHflfliiik 1200 2?!. where r g = GM/c 2 
is the gravitational radius of the black hole with mass 
M. For the mildly relativistic speeds (T ss 1 - 2) consid¬ 
ered here, and jet distances of order of the orbital radius, 
the accretion disk can be well approximated as a point 
source of radiation. The accretion-disk radiation spec¬ 
trum may not be well represented by a blackbody or a 


Shakura-Sunyaev spectrum, as accretion disks display a 
wide range of spectra. This does not affect eq. m , as 
this estimate is based on energy dissipation at various 
radii, which is largely unaffected by disk type. 

Following the procedure in A ppen dix A, again using 
the head-on approximation, eq. ll ATI . for the Compton 
cross section, but now for a point source at the origin, we 
derive the vF v spectrum of Compton-scattered accretion 
disk (CSAD) radiation given by 


f? = 


3e 2 (5f, 


d^t 


1287r 2 a: 2 d 2 


[y- 


de 0 
2 e 


^o(eo) 


, ^e(7) 

d l 




dvie-iy 

where e* = 7 Co<5d( 1 ~ A 4 ) an d 


+ (-T ± -) 2 ] , 

o D 7 e* 2 / 


( 6 ) 


7 min ~ 2 So t 1 



ee 0 (l - y) 


Appendix B gives approximate expressions for the CSSR 
and CSAD spectra in the Thomson and KN regimes that 
reduce the number of numerical integrations, though at 
the expense of accuracy. These simpler expressions have, 
however, the virtue of allowing the beaming factors of the 
various processes to be simply derived. 


3. RESULTS 

We calculate the CSSR spectrum using the standard 
parameters listed in Table 1. The standard jet height x 
is taken to be equal to the mean orbital separation d = 
2.5 x 10 12 cm of LS 5039. We also employ a broken power- 
law distribution for the nonthermal electrons, given by 


N 'e{l) = K [l 1 P 7 7o> 7i) + 7 p #( 7 ; 7 i, 7 2 )] . 

( ? ) 

Normalizing to the total comoving nonthermal electron 
energy 

/ oo 

dj 7 K( 7) 


gives 


K = 


K 

m e c 2 


q—p/ 2—o 2—g\ 2— p 2—p 

r7i (7i — 7o ) . 7i ~ 1 -1 


2 -q 


P~ 2 


• ( 8 ) 


The maximum electron L orentz factor is calculated 
from the well-known dGuilbert et al.l 119831 radiation- 
reaction limit obtained by equating the electron energy- 
loss timescale with the gyration timescale t g = m e c^ / eB , 
which holds provided that synchrotron losses dominate 
Compton-energy losses. The situation is more compli¬ 
cated here, where Compton-losses from the stellar ra- 
diation field are large and KN effects are important 
llAharonian et alJl2005 a':L In Appendix C, an improved 
treatment of this limit for relativistic jets is performed, 
and we quantify the regime where the synchrotron 
radiation-reaction limit hold s. We take 72 = 7 max here, 
with 7 m ax given by eq. ( 1 C 111 and 77 = 1 , giving the most 
optimistic maximum synchrotron frequency. 

Fig. 2 show calculations of the vF v SED of the CSSR 
process. The dependence of the SED on phase of the 
binary orbit is shown in Fig. 2a, and the dependence of 
the SED on various parameters of the binary system at 
<j) = 7 t /2 is shown in Fig. 2b (we now simplify the notation 
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Fig. 2.— The uF u Compton-scattered stellar radiation (CSSR) 
spectral energy distribution (SED) from a microquasar jet when 
nonthermal jet electrons Compton-scatter photons of the high-mass 
companion star. Standard parameters of the system, appropriate 
to LS 5039, are given in Table 1. In particular, the jet height 
above the orbital plane of the system is assigned a standard value 
of x = 2.5 X 10 12 cm, the jet plasma outflow Lorentz factor T = 2, 
and the number index of the electron energy distribution is p = 3. 
(a) Dependence of the SED on phase 0, with 0 = 0 when the 
companion star is nearest the observer (Fig. la). The average 
SED integrated over all phases is shown by the solid curve, (b) 
Dependence of the CSSR SED on x and p with 0 = 7r/2. The solid 
curve shows the SED for the standard parameters, and the other 
curves are labeled with the parameter that has changed from the 
standard parameter values. Effects of 77 attenuation on the SEDs 
are not included here. 


by letting the stellar azimuth angle </>* —></>). In both fig¬ 
ures, the steep lower energy behavior is a consequence of 
the low-energy cutoff Lorentz factor 70 = 100 assumed in 
the electron distribution. The energy of this low-energy 
cutoff, and the cutoff associated with the change in slope 
at 71 = 10 3 , occurs for the blackbody stellar spectrum 
around photon energy 

E t Si 2.7fcBT*7 I 2 <5 D r(l-/?/b t )(l —cos?//) , i = 0,1 , (9) 

provided that the scattering takes place in the Thomson 
regime. (This result is found from the ^-function approx¬ 
imations for Thomson scatering given in Appendix B.) 
Except for the phase-dependent factor (1 —cos ?//), the 
terms in eq. @ are uniquely determined. For the stan¬ 
dard parameters of Table 1, f3 = 0.866, i5d = 2.32 and 
/I* = 0.707, so that the breaks due to Thomson scatter¬ 
ing occur at = 167/(1 — cos tp') eV, in accord with 
the results shown. The breaks in the photon spectrum 
in Fig. 2a occur at lower energies when </ = 7r than when 
<j> = 0 because the 4> = 7r case involves more nearly tail- 
on scattering events that scatter the target photons to 
lower energies than for the more nearly head-on scatter¬ 
ing events with <p ss 0 (see Fig. la). 

The spectral index of the Compton-scattered stellar 
radiation follows the well-known behavior 


in the Thomson regime, as shown in Appendix B, where 
a v is the vF v spectral index 4 . For electron index p = 2 
assumed in the lower branch of the electron spectrum 
in Fig. 2a, the spectrum rises with slope 1/2 in a vF„ 
plot. For the upper branch of the electron spectrum with 
p = 3, the spectrum is flat, as seen for the <f> = t: curve in 
Fig. 2a. For the more nearly head-on scattering events 


Fig. 3. — The vF u CSSR SED as a function of jet height x 
at phases 0 = 0 and 7r, for standard parameters given in Table 
1. (a) Dependence of the SED on x at phases 0 = 0 and 0 = 7r 
shown in the upper and lower curves, respectively, of each pair, (b) 
The value of vF u at 100 MeV at 0 = 0 and 0 = 7r. Effects of 77 
attenuation on the SEDs are not included here, but is unimportant 
below ^ 50 GeV. 


with <j> ss 0, the transition to the KN regime occurs at 
lower photon energies, so that scattering in the Thomson 
regime is never fully achieved. Indeed, the CSSR spec¬ 
tral index at high-energies asymptotically approaches the 
steep KN index a v ~ 1 — p (eq. M)- 

Transition to the KN regime of scattering, as derived 
in Appendix B, takes place over a broad range of photon 
energies centered around 5 

^ m e c 2 5 y> 1 ^ 20GeV 

KN 4x2.70 r(l — /?£,) (1 — cos tjj') (1 — cos?//) 

0 °) 

The three terms making up the middle expression rep¬ 
resent, from left to right, the transition energy for a 
monochromatic radiation field with mean photon energy 
~ 2.7O0m e c 2 , a factor for the dependence on jet speed, 
stellar location and observer direction, and an angular 
factor accounting for the stellar azimuth. The <f> = 0 
and (f) = 7T cases in Fig. 2a have, roughly, cos ip' ~ —0.9 
and cos ip’ ~ +0.8, which accounts for the large ranges 
in photon energies corresponding to the transition to the 
KN behavior. KN effects are quite substantial, however, 
even at lower energies than given by eq. m, because 
of the gradual change in the Compton cross section over 
the range where recoil effects start to become important. 

Because the transition to the KN regime occurs at 
lower values of scattered photon energy at <p = 0 where 
the flux is higher than for the <p = ir case where the 
flux is lower, the CSSR process will produce a softer 
spectra with increasing flux at multi-GeV - TeV ener¬ 
gies. This is the same behavior as inferred from the ef- 
fects of qy absorption on th e emitted radiation spectrum 
lIBottcher fc Dermeil 12005(1 . and so would enhance this 
behavior. 

Besides the rapid transition to the steep KN behavior 
at energies < 100 GeV, Fig. 2a shows that the CSSR 
spectra are strongly modulated in flux as a function of 
phase of the binary system. Also shown in Fig. 2a by 
the solid curve is the phase-averaged CSSR SED, which 
is approximately equal to the <f> = 7r/2 curve (actually 
closer to </> = 7/16). Fig. 2b shows calculations of CSSR 
at cj) = 7t/2 f° r different values of p and x. Due to 


4 At lower energies, a v = 2, corresponding to the low-energy 
emissivit2_S£ectruin2(_2ie , _2fa; mono-energetic electron distribu¬ 
tion felumenthal fc (lonldl l 197(1) . 


° The factor 1/4 is introduced because the domain of Compton 
scatteri ng is defined 1 i v the value of the quantity 47c' rather than 
ye'; see iBlumenthalT'Coiildl J1970D. 
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Fig. 4. — Dependence of the CSSR SED of a two-sided jet on 
the bulk Lorentz factor T of the plasma. The SEDs for the forward 
(f) and reverse (r) jets are shown for V = 1.2 and T = 2 at phases 
(f> = 0, ir/2 , and it (see Tabic 1 for parameters). Only the <j> = 7 t/2 
curve is shown for the T = 2 forward jet; the other phase-dependent 
spectra for this case are shown in Fig. 2a. 


the diminution of the target photon density, the abso¬ 
lute value of the flux decreases rapidly with jet height 
x, as expected. The harder electron spectrum causes the 
peak of the CSSR SED to shift to higher energies, but 
the spectrum still falls rapidly when KN effects become 
important. 

Fig. 3 shows in more detail how the CSSR SED de¬ 
pends on jet height x at phases <j) = 0 and 7 r. For these 
parameters, the vF v spectrum peaks between 100 MeV 
and 1 GeV, as shown in Fig. 3a. If 7 -ray emission from 
microquasars is due to this process, then strong phase- 
dependent modulation is expected if the jet electrons ra¬ 
diate on a distance scale comparable to the orbital ra¬ 
dius. As seen in Fig. 3b, the CSSR flux nulls at (j> = n 
when the inclination angle 9 = arctan (d/x), which takes 
place at x = 5.4 x 10 12 cm for our parameters. This 
happens because the target photons are directed almost 
exactly “tail-on” past those electrons which would scat¬ 
ter photons into our observing direction. This effect will 
be somewhat ameliorated in calculations where the star 
has a finite extent. 

Unlike the 77 attenuation process, which ceases to be 
effective at producing even 10 % orbit al mo dul ation of the 
VHE emission when x ~ 10 13 cm ( (Bottcher fc Dermer! 
nni), the CSSR process will produce significant phase- 
dependent modulation of the ~ 100 MeV - GeV emission 
to much larger distances. The 100 MeV fluxes at = 0 
and cf> = n differ by a factor of 2 at x « 5 x 10 13 cm, and 
by a factor of 10% at x ~ 5 x 10 14 cm. Therefore even 
if the 7 -ray emission is not produced in the inner jet, 
GLAST could still detect significant orbital modulation 
of the 7 rays. The large range of distances that pro- 
du ce si gnificant orbital modulation also makes the claim 
(lMassill2004l) of periodic variability of the EGRET emis¬ 
sion from 3EG J0241+6103 associated with LSI+61°303 
more believable, and suggests that the EGRET data for 
3EG J1824-1514 should be reanalyzed to establish strong 
limits to modulation at the orbital period of LS 5039. 

Fig. 4 illustrates the dependence of the fluxes of a two- 
sided jet for model jets with T = 1.2 and 2. The for¬ 
ward (f) jets are inclined at 9 = 25° to the observer, and 


the reverse (r) jets have 9 = 155°. Even at mildly rel¬ 
ativistic Lorentz factors, the Doppler boosting can lead 
to differences of several orders of magnitude between the 
fluxes from the forward and reverse jets. As shown in 
Appendix B, the beaming factors for Thomson-scattered 
radiation for external photons that enter the jet from 
behind and in front of the jet are oc <5p +p (l — p)( p+1 ^ 2 

and oc 6 q + p (1 + respectively. For the standard 

parameters, p = 0.9063, p = 3, and (3 = 0.866 for F = 2. 
Thus the ratio of scattered fluxes of the forward and re¬ 
verse jets range from « 800 for target photons entering 
from behind, to « 10 ® for target photons entering from 
in front of the jet. Intermediate values are found for stel¬ 
lar photons that enter at shallow angles, as calculated in 
Fig. 4. This shows that CSSR from the reverse jet can 
generally be neglected for even mildly relativistic jets. 

4. COMPTON-SCATTERING LEPTONIC JET MODEL FOR 

LS 5039 

In this section, we apply the proceeding results to LS 
5039 data, showing that 

• It is very difficult to fit the HESS with a CSSR 
model because of the strong curvature produced 
by KN effects; 

• An improbable leptonic model with electrons accel¬ 
erated with maximum efficiency with a distribution 
displaying no spectral breaks is required to fit both 
EGRET and HESS data, assuming that these data 
are representative of the simultaneously measured 
spectrum of LS 5039; 

• The EGRET and HESS data can be separately fit 
with emission that is primarily from the CSSR and 
SSC processes, respectively, which implies signifi¬ 
cant aperiodic variability of the 7 rays. 

We furthermore consider other possible origins for the 
broadband emission from LS 5039, including hadronic 
and combined hadronic and leptonic models, and the 
possibility that the EGRET emission associated with LS 
5039 originates from another source in the LS 5039 field, 
or is contaminated by, e.g., diffuse radi ation . 

Fig. 5 shows non-simultaneous data llAharonian et all 
lama) taken at radio, optical, X-ray, and 7 -ray 
wavelengths, wit h EGRET data points taken from 
[Bosch-Ramon fc Pareded (l2004al) . The HESS data are 
described by a power-law spectrum between 250 GeV a nd 
4 TeV with photon index = 2.1 ± 0.15 llAharonian et alJ 
f2005tb We find that fitting the HESS spectrum with a 
pure CSSR model faces severe difficulties. This is due 
to the strong KN decline and marked convexity of the 
CSSR radiation spectrum at TeV energies. Extremely 
hard (p « 2 ) electron spectra extending to very high 
energies are required for such a model to work, but in 
all cases the CSSR SEDs display strong curvature that 
should be measurable with more sensitive HESS data. 

Assume now that the EGRET and HESS data repre¬ 
sent the 7 -ray spectrum of LS 5039 that would be mea¬ 
sured at a single epoch. In order to fit both the EGRET 
and HESS data simultaneously, we find that a CSSR fit 
to the EGRET data will produce too soft of a spectrum 
in the HESS range. This is likewise the case for an SSC 
model fit to the EGRET data. A leptonic SSC model 
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Fig. 5. — Leptonic jet model fit for the LS 5039 data, using Model 
A parameters given in Table 1. The low and high X-ray spectra 
are from XMM and RXTE observations, respectively. The optical 
data have not been dereddened. This figure illustrates the diffi¬ 
culty that leptonic jet models have to fit the joint nonsimultaneous 
EGRET and HESS data if these data are assumed to represent the 
broadband LS 5039 7-ray spectrum measured at a single epoch. 
Inset shows the fit to the EGRET and HESS data in more detail, 
with the dashed line indicating the small range of power-law spec¬ 
tra that connect the EGRET and HESS data. Total spectrum is 
shown by the solid curve, with components as labeled. 


can fit the HESS data alone, as we show below. But an 
SSC model fit to the HESS data means that the EGRET 
data require a separate explanation. One possibility con¬ 
sidered in Fig. 5 is that the high-energy extension of 
the synchrotron spectrum could produce emission in the 
EGRET band. For the hard synchrotron spectrum re¬ 
quired in this case, it is then necessary to fine-tune the 
joint CSSR and SSC spectra to obtain an acceptable fit 
to the joint EGRET and HESS data (see inset of Fig. 5). 
Model A parameters given in Table 1 are used in this fit, 
noting that the mean magnetic field B < 10 G in order 
to have sufficiently energetic electr ons to scatter target 
photons to multi-TeV energies (eq. eh). 

It is interesting to note that the radio and low-state 
X-ray emission from LS 5039 can be connected by a 
single power-law. This suggests that the radio through 
X/ 7 -ray emission could originate from the same (non- 
thermal jet synchrotron) process, as has been argued to 
be the case for mod els of the LMX Bs XTE J1118+480 
llMarkoff et al.l207nt) and GX 339-4 llMarkoff et all2005) 
on the basis of radio, IR, and X-ray spectral correlations. 
In those models, a broken power-law jet synchrotron ra¬ 
diation spectrum resulting from adiabatic and radiative 
cooling is used to fit the data. A similar type of model, as 
shown in Fig. 5, could in principle fit the radio/low-state 
X-ray data from LS 5039 if the synchrotron cooling fre¬ 
quency lies either below radio frequencies or above ~ 100 
MeV 7 -ray energies. 

Nevertheless, we do not consider this a viable expla¬ 
nation, as it requires acceleration of electrons with max¬ 
imal efficiency (eq. ED) with no strong evidence for a 
synchrotron cooling break from the radio through ~ 100 
MeV 7 -ray regime. It also requires fine-tuning of the 
CSSR and SSC spectra to fit the HESS data. Moreover, 



Fig. 6 .— Leptonic jet model fit for the LS 5039 data, using 
Model B parameters given in Table 1. The heavy and light solid 
curves in the main panel and inset gives the total spectra for Model 
B parameters with x = 10 13 cm and x = 10 14 cm, respectively. 


in this jet model, as well as in the models by Markoff 
et al., the radio emission begins to be self-absorbed be¬ 
low ~ 10 GHz (see Fig. 5) for a comoving jet radius 
rb = 2 x 10 11 cm (Table 1) needed to give the requisite 
SSC flux. The radio data of LS 5039, by contrast, shows 
no evidence for self-absorption. 

Thus if the EGRET and HESS data are to be jointly 
fit, it seems that a leptonic jet model does not work. 
Yet a simple hadronic model, where the emission is due 
to a power-law spectrum of cosmic rays accelerated by 
the outer parts of the jet, which then interact with the 
ambient medium or the d ense matter field of the stellar 
wind itRomero et a,l.ll2003tl to produce secondary nuclear 
pion-decay 7 rays, seems no more likely. The dot-dashed 
line connecting the EGRET and HESS data in the inset 
to Fig. 5 suggests that a single power law with number 
index ~ 2.4 nearly connects the two data sets. But a 
power-law cosmic-ray spectrum produces a convex 7 -ray 
spectrum at a few hundred M eV due to the p ion produc¬ 
tion and decay processes dRomero et all2005l) that would 
miss the EGRET data points which are most significant. 

Any leptonic jet model would predict stochastic vari¬ 
ability, so clearly the assumption that needs to be relaxed 
is that the EGRET and HESS data, taken years apart, 
are typical of a contemporaneously measured data set. 
Leptonic jet models must be highly variable in view of 
the short radiative cooling times of nonthermal jet elec¬ 
trons, which require active mechanisms for accelerating 
the electrons, such as internal shocks from variable rela¬ 
tivistic plasma outflows that form the microquasar jets. 

Fig. 6 shows an example of a fit where the EGRET 
data are due primarily to the CSSR process, and the 
HESS data are due primarily to the SSC process. Here 
we take = 10 11 cm. The spectral components, and the 
heavy solid and dotted curves for the CSSR process and 
the total spectrum, respectively, refer to the case when 
x = 10 13 cm. The light solid curve for the total spectrum 
and the light dotted curve for the CSSR process refers to 
the case when x = 10 14 cm. All other parameters remain 
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the same (we let cj> = 7 r /2 in our fits). The difference be¬ 
tween the two spectral fits are then due to the effects on 
the CSSR component arising from the different distances 
of the jet from the star. We do not fit the radio emis¬ 
sion, whi ch is p roduced scales as large as ~ 1000 AU 

Also shown for completeness is a model accretion-disk 
spectrum with emission that is comparable to the low- 
state X-ray flux, along with the associated CSAD radi¬ 
ation. When x 10 12 cm, the CSAD process makes 
no significant contribution to the 7 -ray emission. It is 
surprising, however, how weak and soft (here we use an 
effective temperature of 0.3 keV) the accretion-disk radi¬ 
ation field has to be in comparison with the flux level of 
the jet radiation, and even here the model overproduces 
the low-state X-ray observations. The use of more real¬ 
istic disk models with hard X-ray emissions could better 
fit the low-state X-ray data, as observed with XMM- 
Newton, as well as X-ray states measured with RXTE 
JBosch-Ramon et aLll2005|) . though these authors favor 
a jet origin for the X-rays in view of the spectral in¬ 
dex/flux correlation and the smooth transition that is 
indicated between the X-ray and 7 -ray components (the 
X-ray high state shown in Figs. 5 and 6 represent the 
historical 1998 high fluxes measured with RX TE, but 
may be contaminated by diffuse emission; see lR.ibo et all 
Il999t IBosch-Ramon et aLll20 05f. 

Fig. 6 shows that if the CSSR process accounts for 
the origin of the 7 rays measured with EGRET, then 
it would exceed the flux measured by HESS at pho¬ 
ton energies of several hundred GeV. This can be re¬ 
versed to make a prediction that when the emission from 
LS 5039 is luminous in the ~ 100 MeV - GeV band, 
then the HE SS spectrum wil l be brighter and softer than 
reported by lAharonian et al.1 ( 1200 5*1. Besides predict¬ 
ing stochastic variability, our model would also predict 
weak or absent periodic variability of the HESS emis¬ 
sion on the orbital timescale associated with the SSC 
emission process. This is because SSC emission does 
not exhibit (for fully isotropized electron distributions) 
the periodic emission signatures of the CSSR process. 
There could, however, be variability of the HESS emis¬ 
sion on the orbital timescale due to 77 effects if the SSC 
emission originated from the inner jet, which would dis¬ 
play the d istinctive featur es that we recently calculated 
(iBottcher fc Dermeill2005tf . Aperiodic variability of the 
7 ray emission could arise from a region as large as ~ 10 
AU, corresponding to the size scale of resolved, radio 
emission from LS 5039. The size scale of the 7 -ray emis¬ 
sion would be difficult to constrain unless sub-hour vari¬ 
ability were measured with HESS or CLAST. Note that 
the < V resolution of HESS corresponds, at 2.5 kpc, to 
< 100 AU. EGRET localized sources in the Milky Way’s 
disk to « 1° at best, whereas GLAST will achieve a lo¬ 
calization by a factor « 10 better than EGRET. 

The situation is, of course, more complicated, and we 
can consider a number of other possibilities to explain 
the broadband spectrum of LS 5039. For example, 

1 . A model involving hadronic cascades from an in- 
ne r jet source might fit t he data. The model 
of lAharonian et al l ( 2005a ) would predict strong 
phase-dependent modulation of the 7 -rays; if this 
is not detected with HESS, this model will probably 


be ruled out. A more probable hadronic model, in 
accord with indications of variability, involves the 
interaction of cos mic-ray accelerated p articles with 
the stellar wind dRomero et al.112003 0. though the 
energetics could be demanding if the cosmic rays 
diffuse rapidly away from the source. 

2. A combined leptonic and hadronic model could ex¬ 
plain the observations. In such a model, the ~ 100 
MeV - GeV 7 rays are due to CSSR with a soft 
(cooling) electron spectrum with p > 3, and the 
VHE 7 rays would be formed by a hadronic emis¬ 
sion component due to secondary nuclear produc¬ 
tion with particles from the surroundin g mediu m 
or st ellar wind, as in the model of IRomero et alJ 
( 12005H . In this case, we would expect variable 
modulated emission at GLAST energies, and non¬ 
varying or slowly varying radiation at TeV energies 
with no periodic modulation. One would also ex¬ 
pect extended VHE emission around LS 5039 as 
the accelerated cosmic rays diffuse away from the 
source. In contrast, LS 5039 emissio n is reported to 
be co nsistent with a point source lAharonian et all 
2005 1; even if the extended VHE emission were too 
weak to detect with HESS, an associated radio halo 
from pion-decay secondaries would be expected if a 
dense target, su ch as a molecular cloud, was nearby 
(fBosch-Ramon et al1ll2005dl . 

3. The EGRET emission could originate from another 
source in the field of LS 5039. Both PSR B1822- 
14 and SNR G 16.8-1 .1, the association proposed 
bv lTorre s et al.l (1200311 . are withi n the EGRET er¬ 
ror box (lAharonian et al . 2 QQII- It is doubtful 
that a hadronic SNR component could produce the 
EGRET emission without also being a significant 
HESS source, so that emission from a pulsar seems 
more likely. Cosmic ray irradiation o f clumped cold 
dust n ot recognized in CO surveys (IGrenier et alJ 
2005lJ could also masquerade as a 7 -ray source in 
the LS 5039 error box. Because EGRET source 
detection depends sensitively on the cosmic-ray 
induced diffuse background, incorrect background 
subtraction for EGRET, which has a much larger 
spatial resolution than HESS, could give an incor¬ 
rect flux for the LS 5039 source. Insofar as LS 5039 
is a VHE source, we would still expect it to radiate 
in the ~ 100 MeV - GeV band, though possibly at 
a flux level lower than reported by EGRET. 

In spite of these possibilities, we think that emis¬ 
sion from a leptonic jet is the most likely explana¬ 
tion for the origin of the 7 -ray emission from LS 
5039 and other galactic microquasars. We agree with 
IBosch-Ramon fc Paredesi (2004a) that the EGRET 7 
rays are likely dominated by CSSR; indeed, a model 
where the EGRET emission is primarily fit with an 
SSC component giv es a poor fit to the HESS data 
dParedes et alJEoOa) . More likely, the VHE emission ob¬ 
served with HESS is mainly SSC radiation. 

5. SUMMARY AND CONCLUSIONS 

In this work, we developed a leptonic model for the 
VHE emission observed from LS 5039 with HESS, where 
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the Compton-scattered flux depends on observer an¬ 
gle 9, bulk Lorentz factor F, and the parameters of 
the binary microquasar system. If the 7 -ray flux of 
a galactic microquasar is due to stellar photons that 
are Compton-scattered by jet electrons, then the 7 -ray 
spectrum will exhibit variability correlated with orbital 
phase of the companion star. The fractional modula¬ 
tion depends strongly on the location of the jet. Phase- 
dependent modulation of the 7 rays can result both from 
the CSS R pr ocess a nd 77 pair-pr oduction absorption 
llBottcher fc Dermeil2005tiDubusl2005l) of jet 7 rays that 
interact with photons of the high-mass star. Periodic 
modulation of the VHE emission with signatures charac¬ 
teristic of the 77 attenuation would provide strong ev¬ 
idence in favor of an inner jet model for the origin of 
VHE 7 radiation from LS 5039, whereas variability of 
the ~ 100 MeV - GeV emission on the orbital timescale 
allows the 7 rays to originate from either the inner or 
extended (d < x < lOOci) jet. 

Our study of the CSSR process for high-mass micro¬ 
quasars yielded a number of results: 

1. The high flux states associated with the orbital 
modulation of the binary microquasar system ex¬ 
hibit stronger spectral softening than the low flux 
states as a result stronger KN effects near orbital 
phase <j> = 0 than at phase <f> = 7 r; 

2. An expression was derived, eq. ©, relating breaks 
in the CSSR spectra for scattering in the Thom¬ 
son regime to breaks in the electron distribution. 
These breaks can potentially be inferred from the 
synchrotron radiation spectrum, and used to test a 
CSSR model. 

3. Another expression, eq. m, was derived that 
gives the photon energy for the onset of KN ef¬ 
fects in CSSR spectra. The KN softening can ap¬ 
pear at energies less than ~ m e c 2 /(4 x 2.70) ~ 
25 GeV/T*(eV), even for mildly relativistic out¬ 
flows, due to the gradual onset of the decline in 
the Compton cross section. 

4. Beaming factors for the CSSR process were derived 
in Appendix B, and it was shown that the principal 
dependence of the CSSR vF v flux f e oc 6^~ p (eqs. 
m and |B9| 1 for an electron spectrum with index 
P- 

5. Ratios of CSSR fluxes from forward and reverse 
jets for a microquasar system were obtained. Even 
for mildly relativistic outflows (r ~ 1.5) of mi¬ 
croquasar jets, the ratios of the CSSR fluxes from 
the forward and reverse jets can be many orders of 
magnitude. 

6 . Allowed ranges of parameters related to maximum 
electron energies permitted by competing Comp¬ 
ton losses and measured ratios of Compton-to- 
synchrotron fluxes were derived in Appendices C 
and D, respectively, showing that a leptonic inner 
jet must be, at least, mildly relativistic. 

7. Significant modulation of the CSSR 7 rays with 
the orbital period of the star can take place for 
emission produced in the extended jet, unlike the 


modulation produced by 77 attenuation, which is 
only important for emission produced in the inner 
jet. 

From these results, we tried to construct a purely lep¬ 
tonic model for the multiwavelength SED of LS 5039. 
The HESS data cannot be fit with the CSSR process 
alone due to the strong curvature of the scattered spec¬ 
tra. If the CSSR process is used to fit the EGRET data, 
then it overproduces the emission in the HESS band. 
Thus, it is difficult for the CSSR and SSC processes to 
fit both the HESS and EGRET data simultaneously. To 
circumvent this difficulty, we considered a model, Fig. 5, 
where the EGRET emission is due to nonthermal syn¬ 
chrotron radiation emitted by an electron distribution 
with a very hard spectrum (p < 2.5) that extends to 
the radiation-reaction limited maximum electron energy. 
Not only is this electron distribution unrealistic by dis¬ 
playing no cooling break, but the model requires fine- 
tuning of the CSSR and SSC processes to fit the HESS 
data. 

Because leptonic jet models should exhibit marked 
aperiodic variability, it is doubtful that the HESS and 
EGRET data are representative of the 7 -ray spectrum 
at a fixed epoch. Abandoning this assumption makes 
it feasible to fit the EGRET data with CSSR emission 
(Fig. 6 ) and the HESS d ata with SSC emission (as in 
the model of lAtovan fc Aharonianl ,ii 999tb This leptonic 
model predicts stochastic variability of the HESS data, 
and enhanced emission at several hundred GeV energies 
when the ~ 100 MeV - GeV emission is in a high-flux 
state. The level of modulation of the AGILE/GLAST 
and HESS 7 -ray flux on the stellar orbital period will 
directly imply the location of the 7 -ray emission site; 
in the former case from the inner or extended jet, and 
for the latter case, from within the inner jet. Absence 
of periodicity of the 7 -ray emission measured with both 
GLAST and HESS means that the 7 rays originate from 
the outer (x > lOOd) jet. Confirmation of a leptonic jet 
model can be made by comparing contemporaneous X- 
ray and 'Y-ray o bservations wit h predictions of correlated 
variability llGunta et al.ll2005lh 

If the jet flow is nonrelativistic, then the jet emission 
must be produced in the extended or outer jet in order 
that electrons can be accelerated to sufficiently high ener¬ 
gies to emit multi-TeV 7 rays (see Appendices B and C). 
Thus, if the outflows forming galactic microquasar jets 
are nonrelativistic throughout the full extent of the jet, 
only weak periodic modulation from Compton scattering, 
and essentially no periodic modulation from 77 absorp¬ 
tion effects, is possible. A mildly relativistic (F > 2) 
leptonic jet model in the inner regions of the jet is com¬ 
patible with the acceleration of sufficiently energetic elec¬ 
trons to make VHE 7 rays throughout the HESS band. 
If HESS data show variability on the orbital timescale, 
it means that nhcroquasars eject relativistic flows, which 
must then be decelerated to nonrelativistic speeds by the 
time the jets are detected with radio telescopes. Jet de¬ 
celeration can be accomplished through radiative drag on 
the stellar radiation field —which makes the 7 rays—as 
well as via interactions with the surrounding medium to 
form shocks that could accelerate cosmic rays. In this 
case, the radiative opacity should be large and a light 
pair jet might be required. 
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The most serious limitation of our study is the use 
of fixed electron energy distributions. In more realistic 
models, power-law distributions of nonthermal electrons 
are accelerated through shock processes and are injected 
and ev olve throu gh ra diative, adiaba tic, and casca de pro¬ 
cesses ijAtovan Hr. Aharoniaul Il999t (Bosch-Ramon et al.l 
j[2005 bah while the jet plasma flows away from the cen¬ 
tral star. To treat such a system, electron energy-loss 
rates associated with Compton interactions in both the 
Thomson and KN regimes need to be derived fo r use in an 
equat ion for electron energy evolution (see lGuota et al.l 
(2005 . for an analytic treatment in the Thomson regime). 
Systems where Compton losses dominate synchrotron 
losses can introduce unusual hardenings in the steady- 
state electron distributi on and synchrotron spectrum 
( Dermer fc Atovanl 1200211 . so that it is not clear under 
what conditions the assumption of a power-law electron 
distribution is valid. Future work will address this ques¬ 
tion. 

In summary, we have developed a leptonic jet model 
for galactic microquasars to fit the EGRET and HESS 7 - 
ray data for LS 5039, complementary to the leptonic jet 


models developed bv iBosch-Ramon et al.l ll2005all . This 
model predicts aperiodic variability, and periodic vari¬ 
ability correlated with the orbital motion of the star that 
depends on the location of the 7 -ray emission site and the 
speed of the jet outflow. Joint observations with AGILE, 
GLAST and air Cherenkov telescopes will quickly reveal 
the actual contemporaneous 7 -ray spectra of high-mass 
microquasars, their variability properties, and whether a 
joint CSSR/SSC model, as proposed here, is correct. 
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APPENDIX 

DERIVATION OF CSSR SPECTRUM 


Quantities in the comoving jet frame are denoted by primes 6 . The invariance of w*(e», H*)/e 6 implies that the 
spectral energy density of the stellar radiation field in the comoving fluid frame is 1 ') = (e'/e*) 3 u£ b (e*, f 2 *), 

where f 2 * = (arccos/x*, </>*), 

e* = re'(l + , /x* = (/u' + /?)/ (1 + /?//) , and <j>' = <j>* — 7 r . 

The reverse transformations are 


Eq. Q becomes 


where 


e' = re*(l — /?/x*) , \j! = (/x* — /?)/(I — /?/x») , and </>* = 0' + n . 

0 e 3 6(n' - /x'J- 4>*) 


u bb {e',Cl') = u 


M* = 


P (1 + /3//)[exp(e*/ 0 ) - 1 ] 

p>*~ 0 


(Al) 


= 0 * — 7 r . 


1 - 0H* 

This can also be written in the form of the specific spectral photon density 

el 5 ((m'- ~ $*) 


n ' bb { e 'i &') — 


m e c 2 [exp(e*/ 0 ) — 1 ] 


(A2) 


The nonthermal electrons Compton-scatter the target stellar photons that intercept the jet. The comoving-frame 
emissivity as a function of scattered photon energy e' s , scattered direction Hj,, and location x is given by 

=m e c 3 e' a j de' j> dQ' d-y j) dn' e (l - /3 e cos tp')n' bb {e', H'; x)n' e (j, Sl' e )( ) (A3) 

llDermer et al ! 19971) . In this expression, n' e { 7 , fl'^d'ydfl^ is the differential number of electrons per unit proper volume 
in the comoving frame with Lorentz factors between 7 and 7 + ^7 that are directed into the solid angle element d£l' e 
in the direction and , 0 / is the angle between the incident photon and electron in the comoving frame. From Fig. 

lb, 

COS ij}' = n'n' e + \/l — /Li' 2 yjl — \x ’ 2 cos(0' — <j>' e ) , 


and dac / de' s dfl' s is the differential Compton-scattering cross section. The speed of the electron is 0 e c, and fi e = 
— 7 2 « 1 for the relativistic electrons of interest in this problem. 

Photons are scattered within a cone of half-opening angle « 1/7 by electrons with 7 > 1. For highly relativistic 
electrons, we employ the Compton-scattering cross section in the head-on approximation, where the scattered photons 
travel in the same direction as the scattering electrons. Thus 


da c 

de' s dn’ s 



n' 


No primes are attached to the electron Lorentz factor 7, however, as this quantity is always referred to the comoving frame here. 
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and 


dag 

de' 


r _i 2e s / e s \2-i „/ / 2yei 

o- [y + V I" (-) J di\t s \ —, 

87 c,: 7<72/ v 7<72/ 27 l + 2 ej 


) 


(A4) 


Ufi neJl flfiStlRlnmeuthal ft; Gonldll 97fitlDermer ft Schlickeiserll 99.811 . where ctt is the Thomson cross section, H{x\ a, b) 
is the Heaviside function such that if ( x; a, b) = 1 if a < x < b and H(x\ a, b) = 0 otherwise, and 

£ / 

y = 1-- and e, = 7 e'(l — (3 e cosft) = 7 e'(l — cos ft) , 

7 

where e' s is the scattered photon energy in the comoving frame. The term a gives the photon energy in the proper 
frame of the electro n, an d def ines the regi me of interaction (Thomson regime for <C 1, and KN regime for a 1). 
Substituting eqs. G3 and 1IA4I into eq. 1X1 gives, after solving the 5-functions, the result 


= U * c7r r l e 's 

where e, = 7 e'(l — cos ft), 




7 2 


C* r _^ 2 6 s / \2i 

e' exp(e*/ 0 )-l ^ + V 7 ^ 2 / + 7 G 2 / 


(A5) 


cos ft = ft ft - ftl- ft? ftl - ft? cos ft 
The equations relating the comoving and observer quantities are 7 


7 = 


(1 + z)e 


e ' V-P , , / , 

— , Us = i- — , and 0 = <j> . 


5d 5d s 1 — ( 3 g 
The limits on the e'-integral implied by the Heaviside function are 

e' 


4 = 


and = 


2 e' c 


27(7 — e^) (1 — cos ft) “ 1 — cos ft 

The uF v spectrum resulting from Compton-scattered stellar radiation (CSSR) is given by 

/?* = <J'Ws, &,) = ^ e'Xj\e' s ,n' s ) , (A6) 

where 5 d is the Doppler factor, eq. 0 , is the comoving volume of the radiating region, and d l is the luminosity 
distance to the source. 

For a uniform emitting region filled with an isotropic comoving distribution of electrons, 

v'-U-. rV\ _ t// n 'eh) _ N eft) 

V b n e ft,fl s )-V b — 

and we obtain eq. 0 for the CSSR spectrum. In this expression, we define 

a = T(1 + j3ftft) , b = 7(1 — cos ft) and g = a/0 . 

The integrals are defined as 


where 


POO 

Ii(u) = / dx 

J U 


Ii = Ii(ui) - ft(u 2 ) , i = 1,3, 

\ C( 2 ) ~u, u< 1 , 

[ (1 + u) exp(— u) , u > 1 


exp (a;) — 1 

ftn) is the Riemann zeta function (£(2) = 7 r 2 /6 = 1.6449 ...), and 

ui = ae'ftO and u 2 = ae' u /& . 

The function I 2 is analytic and is given by 


(A7) 


1 — p~ U2 

I 2 = In (---) . 

' 1 — e~ Ul 


The function 


Hu) = r 

J u 


dx 


;(e x — 1 ) 


(A 8 ) 


(A9) 


Fig. Al compares the approximations for I\(u) and Is(u) with numerical integrations. These approximations, which 
introduce at most « 10 - 20 % errors over a narrow range, are used in subsequent calculations. 

7 The following results can be applied to sources at cosmological distances by replacing e with (1 + z)e in the right-hand-sides of the 
subsequent expressions. 
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x 

Fig. A7.— Solid curves are t he i ntegr als Ii (u), l 2 (u), and 1 3 (u) given by eqs. <A7l - lA9l . and dotted curves are the approximations for 
I±(u) and .Z 3 (w) given by eqs. and respectively. Inset shows the comparison on a linear scale. 


COMPARISON OF APPROXIMATIONS FOR RADIATION PROCESSES 


Simpler expressions for the Compton-scattered photon spectrum from a point so urce can be derived using 5- function 
approximations for the Thomson and KN regimes in the head-on approximation llDermer & Schlickeiser i 1993 ). The 
accuracy of these expressions are indicated by comparing with accurate results obtained through numerical integrations 
of eqs. @ and @. 

In the Thomson regime, we approximate the differential cross section by the expression 

= <J T 5{e' s - 7Ci)U(l - ei) , (Bl) 


where the Heaviside function with a single argument is defined by H(x) = 1 for x > 0, and H(x) = 0 otherwise. The 
Heaviside function in eq. m obviously restricts scattering to the Thomson regim e. 

Following the approach leading to eq. 0 , but now using cross section m and approximating the star by a 
monochromatic point source with differential energy density 

u*(e*,II*;r) = u*5(e* — 2.7O0)5(/r* — /I*)5(</>* — </>*) , u * = -—^ 5 - , 

47IT z C 

we obtain the 5-function approximation to the CSSR vF v spectrum in the Thomson regime, given by 


f T* _ r4 ccr T u*(l - /3a«*) 2 (1 - cos^') 2 r 2 


/* 1 * 

Je = S L 


8ttc1 2 l 




5d 


2.7O0r(l — /3/Z*)(l — cosip') 


(B2) 


Here 


7 = 


2.7O0r5 D (l — /?/!*) (1 — cosip') ’ 


cost/)' = n*n' s - \/l - mCV 1 - M?cos(^' - <£*) , 

and, as before, /4 = (/I* — /?)/(1 — /?/!*) and (j,' s = (n — /?)/( 1 — /3/z). The factor 2.70 arises because the mean energy 
of a photon in a blackbody radiation field is [3C(4)/£(3 )]/cbF — 2.70&;bF. 

For a point source radiation field located behind the jet, /I* —> 1 and cos'tjj' —* /j! s . Replacing 2.700 with the 
dimensionless soft photon energy eo, we obtain 


llDermer et al • 1992 ). where 


pT.pt _ ,6 CTt(1~ H) 2 

h D 32 TT 2 x 2 d 2 L 


nl/£(l-n) 


de a L(e 0 ) 7t^(7t) 


5d \j eo(l - /-0 ’ 


(B3) 


This result is written in a form suitable for integration over the spectrum of an isotropic point source with spectral 
luminosity L(e 0 ). 
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The beaming properties of the scattered spectrum are easily derived from eq. I1B2II . If IV'( 7 ) ex 7 p , then 

fj* oc ^ 5+p)/2 [T(l - /?/*,)( 1 - cosVi')]^ +1 )/ 2 (e/eo) (3 - p)/2 • (B4) 

For a point source behind the jet, / 2 * —> 1, and 5 d( 1 — //) in eq. IB4I) . so that 

/ e T ’ pt oc 5 3+p (l - Ai)( p+ 1 )/ 2 (e/e 0 ) (3 " p)/2 . 

For a point source in front of the jet, /2* —> —1, and 5 d(1 + p) in eq. ill Fill , so that 


/ e T ’ ptf oc 5 3+p (l + At )( p+ 1 )/ 2 (e/eo ) (3 - p)/2 , 

which has the same beaming dependence as an external isotropic photon field llDerme 3 119931 . Because these results 
span the possible locations of the stellar point source radiation field, we see that the principal dependence of the 
beaming factor in the Thomson regime is fj* oc ^D +p (e/eo )*' 3_p ^ 2 multiplied by a n angle-dependent factor; the exact 
beaming factor for the CSSR process in the Thomson regime is given by eq. nm 
In the KN regime, we use the (5-function approximation 


^KN 

de' s 


ln(2e 1/2 ei)<5(e; - 7 )H(d ~ 1) 

o€i 


(B5) 


llDermer fc Sclilickeiseil 1 01)71 iLi fc Wanull 20 f)ll) for the cross sect ion. The CSSR uF„ spectrum of point-source emission 
scattered by a jet in the KN regime is given for approximation IB-bll by the expression 


= it ln[2e 1 ( 2 7KN • 2.7O0r(l - ».)(! - cos,?)]* 


327rd 2 (2.7O0) 2 e 2 


H e- 


5d 


2.7O0r(l-/3/i»)(l -cos^) 


where 


1KN = 


5d 


(B 6 ) 


(B7) 


The KN regime vF v spectrum of an isotropically emitting point-source of radiation located behind the jet is given 
by the expression 

fKN _ c6 3o~T7KN^e(7KN) f°° , L(e 0 ) s r 9 1/2 /i _ „\1 

h i 2 87r 2 r 2 r/ 2 e 2 / , , . “ £ ° ,2 ln i 2e WJ • 

iz57T x a L e Ji/e(i_/i) e o 

The principal dependence of the beaming factor in the KN regime for a point source behind (—) and in front of (+) 
of the jet, from eqs. (IB 6 I and CEJ and the results following eq. (lB4l ) , goes as 


/ e KN oc <5p +p ln[eeo(l ± /i)]e 1_p . (B9) 

Because of the slowly varying logarithmic factor, the Doppler dependence of the beaming factor in the KN regime goes 
as / e KN oc (5 q + p , as in the case for a surrounding isotropic external radiation held dCeorga.nononlos et a, 1112001 1 . 

Fig. B1 gives a comparison of the (5-function (dotted lines) approximations and the accurate (solid curves) calcula¬ 
tions of the CSSR and CSAD processes, respectively. The lower and upper branches of the approximations for the 
CSSR and CSAD processes represent the Thomson and KN 5-function approximations, respectively. Here we treat a 
monochromatic point source of radiation that radiates isotropically from behind a one-sided jet with Lorentz factor 
r = 2. The parameters of the stellar radiation held and accretion disk point source are given by Table 1 parame¬ 
ter study values, though here we use only a single power-law electron distribution between electron Lorentz factors 
7 = 10 2 and 10 8 with p = 2.5. The accretion disk luminosity is 10 36 ergs s _1 , and the accretion disk is assumed to emit 
monochromatically at 1 keV. The emission from this source is plotted as a thermal emitter with effective temperature 
Opt = 1.0/(2.70 x 511). Depending on the accuracy desired, the 5-function approximations may be adequate for 
calculating the spectrum of a point-source radiation held, and is simpler when integrating over non-monochromatic 
accretion-disk radiation spectra. 

It is also useful to compare simple 5-function approximations for the synchrotron and SSC processes with more 
accurate calculations. A 5-function approximation for the synchrotron spectrum is 


where e B = R(G)/4.414x 10 13 
for the SSC spectrum is 


jsyn ^ 


67rd? 


■ca T u B ^N' e ('y s ) , 7 S = y 




G, and -B(G) is the magnetic held in the comoving frame. A 5-function approximation 


fSSC ~ r4 C^T r bu' B 
Je D 12tt d\Vl 



de' e'-'N'J 




(BH) 
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Fig. B8.— Spectra of stellar and accretion-disk radiation fields (long-dashed curves), and accurate and approximate calculations 
of the spectral components emitted by a one-sided microquasar jet. The jet spectral components include the synchrotron spectrum, the 
synchrotron self-Compton (SSC) spectrum, the Compton-scattered stellar radiation (CSSR) spectrum, and the Compton-scattered accretion 
disk (CSAD) spectrum. Accurate calculations for the synchrotron spectrum, the SSC spectrum, and the CSSR and CSAD spectra in the 
head-on approximation are shown by the solid curves, whereas 5-function approximations are shown by the dotted lines and curves. The 
parameters of the calculation are given by the parameter study values in Table 1, except that here the electron spectrum is a single power 
law with spectral index p = 2.5 for 10 2 < 7 < 10 8 , and <f) = 0. 


llDermer fc SchlickeiserHl2002 '). where u' B = B 2 /8tt and = 47rr^/3. The only provision made in this approximation 
to account for scattering in the KN regime is to terminate the scat tering when ye' > 1. 

For the broken power-law electron distribution, eq. 0, eq. (IB 111) gives 


where 


fSSC * 3 


AJ rp Mg 

D 167r 2 d 2 r 2 


K 2 (—) 3/2 ['Y? q ~ p) l 1 + 7 r P (X 2 +1 3 ) + 24] , 




r -dL\- q l 2 ] r min ( £ ' 

1 y e B } [ max(e'/ 7l 2 ,e s7o 2 ) J 


(B12) 


2 2 = 4 q —{[ min ( e ' 1 > e '/l\i e B7i] (P ,)/2 -[max(e'/ 7 |, e s7 ^)] ,,/z } 


1 (p-9)/2. 


q f P/ 2 f '-9/2 

I 3 = ^i-{[; 

q-p 


nnn(e / 1 , e'/7o> e Bll] (? p)/2 - [max(e'/7i , £b7i)] (? P)/2 } , and 


z,=(T)-^i„ [ mi ° ( 77'7 ,eB J l) i 

es max(e7 7 J,es 7 f) 

(see also lTavecchio et~ .!; 1998 ). In deriving this expression, we assume that the mean escape time of a synchrotron 
photon from the uniform spherical emitting region is 3rf,/4c. 

The quality of the 5-function approximations for the synchrotron and SSC processes is shown in Fig. Bl. Note 
that synchrotron self-absorption plays a role at the lowest synchrotron photon energies, and is not considered in the 
approximate curve. 


ALLOWED PARAMETERS FOR LEPTONIC JET MODEL FROM RADIATION REACTION LIMIT 

We first generalize the well-known expression for the maximum ele ctron Lorentz factor obtained by comparing the 
optimal Fermi acceleration rate with the synchrotron energy-loss rate dCnilhert, et, a,l.llll983t iRaehen ro 1M 

to include Compton losses on the CSSR field. Besides treating KN effects on the energy-loss rates (lAharonian et alJ 
I2005alh our results also correct for relativistic jet motions on the Doppler-boosted stellar radiation field. 

The electron gyration frequency ui g = eR/m e c 7 in a magnetic field with mean intensity B. Note a number of 
underlying assumptions in this treatment: that the mean magnetic field is uniform throughout the emitting volume, 
and that it is randomly oriented. The acceleration rate Lo aC c in Fermi processes cannot exceed the gyration frequency 
because an electron gains at most a fraction of its energy when executing a single gyration in first-order shock or 
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second-order stochastic Fermi processes. Thus u> aC c = V^'g i with 77 < 1. The fractional synchrotron energy-loss rate 
u> syn = | — p S yn/l\ = <JtB 2 j/& mn e c which, when equated with the acceleration rate, gives 


r {max — V 


1/2 


67re _ 1.2 x lO 8 ^ 1 / 2 


(Cl) 


y/mG) 

and a maximum synchrotron frequency of hv max = S D (heB/m e c ) r ) 2 /(1 + 2 ) = 160t7<$d/(1 + 2 ) MeV. 

Scattering in the Klein-Nishina regime is important when 47 c' = 47[2.7O0T(1 — Pp*)] > 1 (compare eq. |10|h that 
is, when 7 > 5 x 10 4 /[/cB , 7*(eV)r(l — Pp*)\. The fractional energy-loss rate in the extreme KN regime for a blackbody 
radiation field with dimensionless temperature 0 is 

wkn = I—— ln(O.552 7 0) (C2) 

7 16 A£ 7 


( Blumenthal fc Gouldl Il97f1l . where Ac = Ac/27r = h/m e c = 3.86 x 10 11 cm is (also called) the electron Compton 
wavelength. 

The energy density it* of a stellar blackbody radiation field in the comoving jet frame has to be corrected by a 
graybody factor that accounts for ( 1 ) the dilution of the radiation field due to the distance of the jet from the star, and 
(2) the reduction in energy density due to the bulk relativistic motion of the jet. Recalling that the energy density at 
the surface of a blackbody is 1/4 times the energy density in the interior of a blackbody cavity, the relation between 
the energy densities is therefore given by 


4 = j(—) 2 r 2 (i -pp*) 2 u bb (0) 
4 r 


(C3) 


(e.g., iDermer fc Schlickeised 1200 2). where Ubb(Q) = 7r 2 m e c 2 0 4 /15Ac is the energy density of a blackbody stellar 
radiation field with dimensionle ss te mperature 0 . 

The graybody factor from eq. d therefore gives the maximum electron Lorentz factor 7 max by solving the equation 

a TB 2 lm ax , (R*\ 2 ^ 2,1 Or-. ccr T Q 2 ,„rn Cr :o-. Or-. M ( C4 ) 


eB 


V 


+ (—) r 2 (i - PP*) 2 —^r e 2 in[o.552 7roox r©(i - pp*)). 


64Ac 


m e c 67 rm e c ' r 

which holds when 7 max [T0(1 — PpP)]^ 1 . The maximum synchrotron energy in the comoving frame is hi/' m 


(heB/m e c)'y/ nax - Eq. IC14I can be rewritten as 


'Ymax V 


9 B c 


3nctf B cr 2 r i?*r(l — /3/u*)©-, 2 


32 


(^) 2 [ 


ln[O.552 7mQ;c r0(l - pp*)] 


(C5) 


4 ctfB 

where cn/ = 1/137 is the fine structure constant and B cr = m 2 c 2 /eh = 4.41 x 10 13 G is the critical magnetic field. 
Again, this expression holds as long as the argument of the logarithm is much greater than unity. 

Because of the similar dependence of the acceleration rate and the fractional KN energy-loss rate found in eq. 
Compton-scattering in the KN regime will be unimportant to limit electron acceleration when 


JL* 


1 7 tB c 

24^B 


= 0.07 


T*(eV) T( 1 - pp*) 

VW( G) 


(C6) 


here we have taken the square root of the logarithmic term ss 2 . 

A limit on the range of values of B can be obtained by noting that the emission of photons with energy FA, requires 
electrons with Lorentz factors 


7 


> 


FA, 


2 x io 7 f; 


10 


S D m e c 2 S D 

where Bio = B 7 /10 TeV. From Eq. lICHl . 

B( G) < 3877<5 2 /E\ 0 . 

With eq. 1 IC 6 I 1 . this implies an allowed range of magnetic fields, given by 


5X 10 3 T 2 (eV) [ R * T{1 < B{ G) < —2 


38 ? 7<5 


(C7) 

(C8) 

n l ~ ~ e 2 ’ 

V r 12 -^10 

where r ±2 = r/10 12 cm and p * = \J 1 — cl 2 /r 2 . For a nonrelativistic outflow, as treated bv lAharonian et j,i. (12005afl . 
77 —> rjop 2 = (0.2) 2 ?7o(/ 3/0.2) 2 , with 70 ^ 1- The allowed range becomes 


V ■ ■ ■ - r . - F .'| 0 

for a leptonic jet model to apply to the observed VHE emission from a microquasar. 
For the parameters of LS 5039 given in Table 1, 

0.024 [r(l — Pp*)} 2 38 / 7 51 

~ -°( 0 ) ro tt 2 ) 


0.6 


Vorl 2 (P/0.2) 2 


< B{ G) < 


1.4t7o(/3/0.2) : 

E 2 
-‘Ao 


(CIO) 
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If the jet is nonrelativistic, eq. itTTol) circumscribes the mean comoving magnetic field of a leptonic jet model for LS 
5039 to a narrow range of values unless r \2 1 even for the optimistic case of 770 = 1, in which case the orbital 

modulation of the 7 -ray emission due to 77 attenuation and, to a lesser extent, the CSSR process will be small. Even 
a mildly relativistic jet (1? > 1.5) opens up a much larger allowable parameter space to permit the sou rce of emission to 
originate from the inner jet. Note that the inferred jet speed from radio observations ijPa,redes et, a,l.1120021 of LS 5039 
may not be representative of the jet speed near the source due to bulk deceleration caused, for example, by radiative 
drag. 


ALLOWED PARAMETERS FOR LEPTONIC JET MODEL FROM RATIO OF COMPTON TO 

SYNCHROTRON FLUXES 


A further restriction on parameters for a leptonic microquasar jet model can be obtained by noting that the ratio 
of the peak vF u fluxes for the Compton and synchrotron components, denoted by p 1 is related to the ratio of the 
comoving photon it* and magnetic field u' B energy densities by the relation 


P = 


' e pk,C 


' €pk,s 


(Dl) 


(compare lSikoralll997l for blazars). This expression, which neglects angular scattering effects, holds for scattering in 
the Thomson regime and so would best be applied to emission observed in the GLAST band. Taking 


u* = 




47 xr 2 c 


(compare eq. 1C31 : ctsb = n 2 m e c 3 k B /[60%Q(m e c 2 ) 4 ] ), we find 


r > 




using eq. O- For the parameters of LS 5039, 

_ > Q-6^0 r(l-flz.) _ 15^0 

ri2 ~ ryvWlOO % 77o(/?/0.2) 2 v^7IOO ’ 


(D2) 


(D3) 


where the last expression holds for a nonrelativistic outflow. Application of this constraint implies that we properly 
identify the peak nonthermal jet synchrotron flux, which is not necessarily obvious from the broadband data (see Fig. 
5). Choosing a value p « 10 2 suggested by Fig. 6 data, we conclude that emission from the inner jet is possible for a 
mildly relativistic leptonic jet model, but not for a nonrelativistic jet at the base of the microquasar. 
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TABLE D1 

Standard values used in parameter study and model fits to LS 5039 


Quantity 

Parameter 

Study 

Model A 

Model B 

d L (kpc) 

3.0 



T,(K) 

39000 



L*(ergs s - 1 ) 

7 x 10 38 



iJt(cm) 

6.5 x 10 11 



d( cm) 

2.5 x 10 12 



9 (inch) 

25° 



r 

2.0 



Sn 

2.32 



B( G) 

1.0 

0.8 

0.8 

Tj,(cm) 

10 11 

2 x 10 11 

10 11 

VF'(ergs) 

10 38 

2 x 10 37 

2 x 10 38 

V 

3.0 

2.25 

3.0 

<1 
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